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A new class of ID discrete nonlinear Schrodinger Hamiltonians with tunable nonlinerities is intro- 
duced, which includes the integrable Ablowitz-Ladik system as a limit. A new subset of equations, 
which are derived from these Hamiltonians using a generalized definition of Poisson brackets, and 
collectively refered to as the N-AL equation, is studied. The symmetry properties of the equa- 
tion are discussed. These equations are shown to possess propagating localized solutions, having 
the continuous translational symmetry of the one-soliton solution of the Ablowitz-Ladik nonlinear 
Schrodinger equation. The N-AL systems are shown to be suitable to study the combined effect of 
the dynamical imbalance of nonlinearity and dispersion and the Peierls-Nabarro potential, arising 
from the lattice discreteness, on the propagating solitary wave like profiles. A perturbative analysis 
shows that the N-AL systems can have discrete breather solutions, due to the presence of saddle 
center bifurcations in phase portraits. The unstaggered localized states are shown to have positive 
effective mass. On the other hand, large width but small amplitude staggered localized states have 
negative effective mass. The collison dynamics of two colliding solitary wave profiles are studied 
numerically. Notwithstanding colliding solitary wave profiles are seen to exhibit nontrivial nonsoli- 
tonic interactions, certain universal features are observed in the collison dynamics. Future scopes 
of this work and possible applications of the N-AL systems are discussed. 
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I. INTRODUCTION 

As is well known, different nonlinear models can pos- 
sess spatially localized solutions for solitary waves 0, g, ||. 
In many cases, the solitary waves are analyzed in the 
frame work of integrable models which, however, describe 
realistic physical systems with certain approximations Q . 
Two important examples of integrable nonlinear equa- 
tions are the nonlinear Schrodinger (NLS) equation, and 
the sine-Gordon (s-G) equation. While the first one is 
known to form " dynamical solitons" , the last one yields 
"kink" and "antikink" solutions. These are also called 
" topological solitons" |IJ . Dynamical solitons arise from 
the acute balance of nonlinearity and dispersion. The ori- 
gin of topological solitons is the balance of nonlinearity 
and constraints from topological invariants [[i] ||. Along 
side these continuous equations, a pioneering example of 
integrable discrete differential equation is the Ablowitz- 
Ladik equation. This is often referred to as the integrable 
discretization of continuous NLS (ALDNLSE) §, 0, §. 
On the other hand, the standard discretization of NLS 
gives the nonintegrable discrete nonlinear Schrodinger 
equation (DNLSE)[|, |l0|, |ll|. 

The application of these continuous integrable equa- 
tions in physics is quite extensive. The NLS, which 
has both dark and bright solitons, was first used to 
understand self-focussing and self-defocussing of carrier 
waves [p"2|. This equation later made inroads in the field 
of biological energy transport in the guise of " Davydov's 
soliton"|| [13| |l|, Quite recently, the NLS is in- 

vestigated in relation to other phenomena, like Bosc- 
Einstein condensates in standing waves, discrete solitons 
and breathers with dilute Bose-Einstcin condensates, and 



the formation of unstable solitonic complexes ]i~6|, fl7| . 
The mathematically very rich s-G equation has also many 
applications ||. It is used as a model for the propa- 
gation of dislocation in crystals and also for the prop- 
agation of magnetic flux in a long Josephson junction 
transmission linejl], ||). On the other hand, the util- 
ity of the ALDNLSE in the analysis of physical prob- 
lems is rather tortuous. Consider for example, physi- 
cally motivated models, such as coupled nonlinear atomic 
strings with onsite or intersite anharmonic potentials fli"8|], 
array of coupled optical wavegui des ||19| , proton dynam- 
ics in hydrogen-bonded chains^, Q |20|, pl| , the Davy- 
dov and Holstein models (DHM) for transport in bio- 
physical systems |2| [|, ^l), and so on. In these mod- 
els, either the ALDNLSE does not appear at all or it 
does not appear in its pristine form. However, to study 
the soliton dynamics perturbatively in these models, the 
ALDNLSE is the appropriate choice for the zeroth or- 
der approximation jl8| [l9, |2l], ^|. This clearly shows 
that the ALDNLSE is an equation of great deal of sig- 
nificance in nonlinear science. For other applications of 
this equation we mention the following examples. The 
dynamics of low frequency and high frequency intrinsic 
localized modes in nonlinear lattices can be described to 
a good approximation by the ALDNLSE pi, 24 1. Again, 



in the study of dark and bright excitons in systems with 
exchange and dipole-dipole interactions, it is shown that 
in some limiting cases the evolution equations resulting 
from model Hamiltonians are reduced to the exactly in- 
tegrable ALDNLSE |^| . As another example, we cite 
the one-dimensional Frolich model of exciton(vibron)- 
phonon interaction. This model can be approximated 
by the ALDNLSE, provided we are interested in the dy- 
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namics of large width and small amplitude dynamical 
solitons [^l]. One other important use of the ALDNLSE 
lies in the numerical integration of the NLS. This is done 
to avoid any numerical instability problemj(|. 

Any small perturbation may break the integrability 
so strongly that solitary waves can be unstable and in 
the extreme case may altogether disappear. The nonin- 
tegrablity in an otherwise integrable nonlinear equation 
can also give birth to internal modes, often called "shape 
modes" in solitary waves. These internal modes can dras- 
tically modify the soliton dynamics pq|. The physical 
origin of integrability breaking terms can be many. For 
example, this type of terms can arise from taking into 
consideration the effect of thermal fluctuations and the 
possible absence of order in the system^, |27], ^8], p9[ . 
However, the most important one is the discreteness of 
the underlying space. Consider for example, the Frenkel- 
Kontorova (FK) model. It is the spatially discretizd s- 
G equation. While the s-G equation which is the long 
wavelength approximation of the FK model is integrable, 
the FK model in its generic form is nonintegrable. This 
nonintegrablity arises from the discretization reflecting 
the discreteness of the space The other interesting 
as well as important example, of course is the already 
mentioned DNLSE§, [l0|, |j. So, relevant in this con- 
text is the study of an IN-DNLS, and also the modified 
Salerno equation (MSE)§l], ||. The IN-DNLS is a hy- 
brid form of the ALDNLSE and the DNLSE with a "tun- 
able" nonlinearity. On the other hand, in the MSE, the 
usual DNLSE is replaced by a modified version of the 
DNLSE, the ADNLSE which involves acousitc phonons, 
instead of optical phonons in condensed matter physics 
parlance |22|. These equations are studied with the prime 
objective to understand the role of lattice discreteness 
as a mechanism for the collapse of moving self-localized 
states to stable, but pinned localized states^]], |3l|. We 
further note that an IN-DNLS, which can be contin- 
ually switched from the DNLSE to the ALDNLSE by 
varying a single parameter, is also studied to investigate 
the discreteness induced oscillatory instabilities of dark 
solitonsp2[. 

There is at least one continuous nonlinear equation, 
"</> 4 equation", which is nonintegrable, contrary to its 
"cousin" s-G system. This </> 4 equation can have either 
solitary wave solutions or " kink- like" solutions with per- 
manent profile. However, these solutions do not have 
the simple collision properties of solitons. These solu- 
tions can bump, lock or annihilate each other. In addi- 
tion these always emit some oscillatory disturbances or 
radiation in the course of a collision. Furthermore, the 
head-on collisions of kink and antikink pair of solutions 
of <fi 4 equation will settle either to a bound state (bion) 
or to a two-soliton solution. This settling down is found 
to depend fractally on the impact velocity. We further 
note that all solutions </> 4 equation have profound physical 
as well as theoretical significance^, [54], |35|, [36], |37) . In 
continuation, we mention that the appearance of fractal 
structure in the kink-breather interaction is investigated 



using the FK model |MJ. Similarly, a fractal structure in 
solitary-wave collisions for the coupled NLS equations is 
reported. This structure is observed in the separation 
velocity versus collision velocity graph Q. 

On the contary, to the best of my knowledge there is 
no known analogy of (/> 4 type of equation in the discrete 
case. Here I propose an extended nonintegrable version 
of the ALDNLSE, which has a "tunable" nonlinearity in 
the intersite hopping term. At the same time, the form 
of nonlinearity is such that it can allow solitary wave like 
solutions, as seen in </> 4 equation, and its suitable parallel 
in the linear regime is one-dimensional correlated disor- 
dered systems p9|, |40[ |. Inasmuch as this nonlinearity is 
in the intersite hopping term, it serves two important 
purposes. First of all, this extra dispersive correction 
to the ALDNLSE will try to destroy the Ablowitz-Ladik 
(A-L) soliton by dispersion. So, by varying this term 
we can investigate the effect of dispersive imbalance on 
the maintance of the moving solitonic profile. It is rel- 
evant at this point to note that both the IN-DNLS and 
the MSE investigate the competition between the on-site 
trapping and the solitonic motion of the A-L solitons. In 
case of the MSE, it is found that the narrow A-L soli- 
tons, having width smaller than the critical width will 
get pinned by the lattice potential^]] ^2|. Similar re- 



sults are obtained numerically from the IN-DNLS 1 31 1. 
Secondly, since the extra dispersive term in the proposed 
equation breaks the integrability of the ALDNLSE, the 
dynamics of the A-L solitons will not be transparent to 
the lattice discreteness. So, along with the IN-DNLS 
and the MSE, this model also gives an opportunity to 
study further the effect of Peierls-Nabarro (PN) poten- 
tial on the dynamics of solitons fl3l[ El). We further note 
that the DHM yields a very complicated nonlinear dis- 
crete differential equation plj. So, to gain a better un- 
derstanding of the DHM, a somewhat simpler but phys- 
ically relevant system needs to be considered pl|. This 
is another important motivation of this study. It is in 
fact noteworthy in this context that a slightly modi- 
fied ALDNLSE has been studied as a plausible model 
for dynamical self-trapping in discrete lattices [Q. As 
for further motivation of this study, we note that one 
rapidly emerging field in nonlinear dynamics is the study 
of the solitary-wave interaction in nonintegrable nonlin- 
ear models. The emergence of this field is due to the pos- 
sibility of observing many of the predicted effects exper- 
imentally including the soliton energy and momentum 
exchange || |§ || ||, |§ 11 @- The model, 
that is to be proposed here allows with an opportunity 
to study the head-on collision of scalar lattice solitary 
waves. Finally, nonlinear phenomena are quite a common 
occurence in all branches of natural science. But, our 
theoretical understanding of nonlinear equations, partic- 
ularly that of nonlinear discrete differential equations is 
quite limited. So, when viewed from this angle, another 
important dimension gets added to the present study. 

The organisation of the paper is as follows. In what fol- 
lows, I first propose a set of discrete nonintegrable Hamil- 
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tonians, and derive the equation of interest from the req- 
uisite Hamiltonian. I study then analytically the exis- 
tence of moving solitary wave like solution in this equa- 
tion. Since this equation is a perturbated Ablowitz-Ladik 
equation, I study this problem further using a standard 
perturbative method to find trapped and moving soli- 
tons. 1 consider next the formation of stationary local- 
ized states and the effective masses of these states. Sub- 
sequently, I present some numerical results to substanti- 
ate the analysis. I also study the collison of two solitary 
wave profiles numerically using the equation. Findings 



from my study are summarized at the end . In this sec- 
tion I also discuss the applicabilty and future scopes of 
the present work. 



II. A GENERAL DERIVATION OF THE 
ABLOWITZ-LADIK CLASS OF NONLINEAR 
EQUATIONS 

We consider here the following Hamiltonian, H. 



h = 



'M+l 



J n+1 H'ni 



1 2v 



n j 



!(£-J)5>[l + A 



(2.1) 



ft is understood that the Hamiltonian, H in Eq.(2.1) de- 
scribes a conservative system in which a quasiparticlc 
(exciton, vibron etc.) moves in an one-dimensional chain. 
To obtain the time evolution equations for generalized 



coordinates, {(p m , '/'ml; we make use of the nonstan- 
dard Poisson bracket relationship among these general- 
ized coordinates as given in Appendix A. After defining 
<p m = (-l) ro ip m , we get for m = 1,2, ... 



Fj{tPm,g J o) = g J i^m^m+aj + 1p m +*j Ipm) 1pm+aj 

a = ±1 

I 

FWm,X) = (1 + A \ip m \ 2 ) ]T Fjilpmid) (2-2) 

j = 1 

i 1p m - 2J 1p m + J (1 + A IV'ml 2 ) (lpm+1 + V'm-l) = 2 V \lp m \ 2 1p m - J r (l/'m,A). (2.3) 

I 



In T{%p m , A), arguements, g®, j — 1, 2, • ■ ■ , I have been 
suppressed, and J > 0. When I — 1, in Eq.(2.1) 
and consequently in Eq.(2.3) we have just an extra non- 
linear nearest-neighbor coupling in hopping with a cou- 
pling constant, g$. This coupling is assumed to arise due 
to quasiparticle-phonon interaction Similarly, when 
I = 2, we have both nonlinear nearest-neighbor and 
next-nearest-neighbor couplings. Coupling constants are 
<7q and g 2 , respectively. So, in this model coupling con- 
stants are assumed to depend on the distance between 
the sites involved. For any arbitrary I then, a given site, 
over and above the standard linear nearest-neighbor cou- 
pling, is coupled through nonlinear hopping to I consec- 
utive sites in general on both sides of it. Of course, the 
origin of such coupling is assumed to arise due to the 
heuristic generalization of quasiparticle-phonon interac- 
tion beyond the nearest neighbor. We further note that 



in this model also 

oo 

N = ^ ln[l+ A iV'ml 2 ] = Constant (2.4) 

m— — oo 

as in the original Ablowitz -Ladik equation (see also 
appendix A). So, if A \ip m \ 2 -C 1 for all m, we get 
l^ml 2 ~ Constant. We now consider the following 
limits of Eq.(2.3). (1) When g° Q = for all j, we have 
the Salerno Equation in the quantum version. In the 
classical domain, it has been nomenclatured as IN-DNLS 
and the formation of staggered localized states from this 
equation is studied j3l| H|. (2) When v = as 
well gjj = for all j, we have the standard Ablowitz- 
Ladik equation. This equation is known to have a single 
soliton solution^, (3) When A = = v and 
also I — 1, we have the model where exciton-phonon 
interaction affects only the hopping between nearest- 
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neighbors. Of course, to obtain this particular non- 
linear mathematical form of interaction, it is assumed 
that the lattice relaxation is faster than the quasipar- 
ticle dynamics This is a very standard assumption 
which is used in the Davydov's soliton [jl3| and also 
in the Rashba and Toyozawa mechanism for the forma- 
tion of self-trapped exciton|l7j This equation 
is also studied for single soliton solution, using pertur- 
bation method Q. (4) The equation with v = will 
be refered to here as the nonintegrable Ablowitz-Ladik 
(N-AL) equation. Consider again the situation where we 
take I = 1 in the N-AL equation and further ignore al- 
together in Eq.(2.3) the term having g^ A. We have then 
a model equation which describes a truncated version of 
the model, in which cxciton(vibron)-phonon interaction 
affects both site-energy and hopping. We note that a 
perturbative analysis of the full model is done for one 
soliton solution pl[. Here our plan is to analyze various 
aspects of the N-AL equation. 

Before we proceed furthere, we define t = Jt, ^> m — 

V~X "0m, m £ Z, and gj = j&j. We further note that 
- -^^m,9j), j = 1,2, ••• I. This, 
in turn yields ^"(^f, A) = ^ J"(* m ,l), m £ Z. So, 
in the limit of v = 0, we have from Eq.(2.3) and m £ Z 

i * m - 2 * m + (1 + |*m| 2 ) (#m+l + *m-l) 

= -F(9 m ,l) (2.5) 

and Eq.(2.4) also remains valid with A = 1 and ip m 
replaced by \t m . Again, T($! mi l) should have in its 
arguement g\ , <?2 , • • • , gi , which have been suppressed 
for convenience. We note that Eq.(2.5) possesses a re- 
flection symmetry. To understand this, we consider the 
transformation, \I> m — > (—1)™ ^ m exp[— Air]. The 
equation will remain invariant under this transformation 



2(1 + Ql) W m (r) 



For convenience in further discussion, we consider only 
one term, say the 1-th term of the sum in Eq.(3.5). We 
note that for cosfc/ = 0, permissible values of k are 
k = 7T - (2^tl) (f), ji = 0,1,2,---,(2Z-1). So, 
there are 21 permissible values of k for k £ [— tt, tt]. We 
now state the following results. 

(1) ForZ = 1,^(^,1) = for to eZ, if|fc| = f. 
So, for this case Eq.(3.1) with k = ±^ are the solitary 
wave like solutions of Eq.(2.5). We further note that 
solitary waves have the maximum possible speed. 

(2) I is odd as well only odd values of j in the sum in the 
Hamiltonian, H are permissible. Also, in this case, only 



provided r — > — r and gj — > —gj, j € I- Again, if 
jF(^ m , 1) = 0, to 6 Z, the equation becomes self-dual 
under this reflection transformation. 



III. PROPAGATING SOLUTIONS OF EQ.(2.5) 

In order to gain an understanding of the nature of the 
solution of Eq.(2.5), we consider the case where gj = 
for all j. As mentioned earlier, in this limit we have the 
integrable Ablowitz-Ladik discrete nonlinear Schrodinger 
equation(ALDNLSE)§, ||. The exact one-soliton solu- 
tion of the ALDNLSE, m € Z is 

t0 , , sinh fi , 

W m (r) = r~f 7 exp[ifc(m - a;) - zaj 

cosh[/i(TO — x)\ 

= Qm( T ) exp[i/c(m — x) — ia] (3-1) 
with the following equations for the soliton parameters: 

(i = 0, jfc = 0, a = u>, (3.2) 

uj = 2 [1 — cosh /i cos k], (3-3) 
2 

x = — sinh/x sinfc. (3-4) 

So, for each /x there exists a band of velocities [see 
Eq.(3.4)] at which the localized state or the one-soliton 
state can travel without experiencing any PN pinning due 
to the lattice discreteness [pi] El. We also note that for 
a given k, fj, £ [0, oo]. Mathematically, one soliton solu- 
tion of the ALDNLSE describes two parameters, namely 
k and /i, family of curves. So, even if we pin the value of 
one of those parameters to a prescribed value, Eq.(3.1) 
will still be a solution of the ALDNLSE. 

We now consider J 7 ^^, 1). Introducing Eq.(3.1) in 
Eq.(2.5) we get 



I 

|fc| = -| is allowed. So, Eq.(3.1) with these particular 
values of k are the solitary wave like solutions of Eq.(2.5). 

(3) I is even as well only even values of j in the sum in 
Eq.(3.5) are permissible. In this case, there is no permis- 
sible value of k. So, Eq.(2.5) will not have any solitary 
wave like solution. 

(4) / is arbitrary and j takes odd values with at least 
one even value. In this case, even if / is odd, there is no 
permissible value of k. So, Eq.(2.5) has no solitary wave 
like solution. 

(5) There is only one term, say Z-th term in the sum in 
Eq.(3.5). In this situation, Eq.(3.1) describes the solitary 



I 

= X! COSfej [(Qm+j + Qm-j) cosk j + KQm+j - Qm-j) sinfcj]}. (3.5) 

3 = 1 
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wave like solution of Eq.(2.5) with 21 permissible values 
of k which are already given. 

We also note that when \gj\ « 1 for j — 1, 2, ■ • • Z, 
Eq.(2.5) then can be treated as a perturbed ALDNLSE. 
In this situation we can use the standard perturbation 
theory to investigate the soliton dynamics from Eq.(2.5). 
This aspect is considered next. 



G 2 (k,l,fi) = Gi(k, I. fx) 



dG(l,k,n) 
dk 







sinh fi dS(fi,x) 
8 _G(M,/0^ 

sinh/i dcosk 
2 G »(M,A»). 



(4.1) 
(4.2) 

(4.3) 



IV. THE PERTURB ATIVE SOLITON 
DYNAMICS OF EQ.(2.5) 



The method is amply discussed in the literature) 21 
|22| , |5C| . In this method, it is assumed that the zeroth 
order solution is the standard one-soliton solution of the 
ALDNLSE. It is further assumed that the effect of pertur- 
bation on the soliton dynamics can be adequately taken 
into consideration by allowing four parameters, namely 
x, k, fi and a to vary adiabatically in time, t. Applica- 
tion of the method to the N-AL equation gives for fi, k 
and x : 



We first note thar the famous Poisson sum formula is 
used to obtain S((i, x) and G(fc, I, mM^i)' ^ n * ne situation 
where cosfcj = for all j, both G(k,l,(i) and 
are zero. In this case, it is easy to see that fx, k and x are 
given by Eq.(3.2) and Eq.(3.4) respectively. Of course, 
for x only certain values of k for k S [—if, if] are allowed. 
Notwithstanding this, this case is as expectedly similar 
to the one-soliton solution of the ALDNLSE. We consider 
next the problem of an effective Hamiltonian. 



S(jj,,x) 
G(k,l,(i) 
Gi(k,l,(i) 



E 



7V~ * 



. sinh ^ 

S — 1 fl 



cos 2nsx 



Ecos 2 kj 
. = 1 sinh m 

, sinh 4 /i . , „ 

4 ^- kt coth^i - 1 - 2 S(u,x)] 

M 



A. An effective Hamiltonian 

In order to derive an effective Hamiltonian for the dy- 
namical system described by Eqs.(4.2) and (4.3), we mu- 
tiply these two equations by x and k respectively. Now, 
substracting the first one from the second one, we get 



. sinh/i dcosk , sinh 4 /! . , „ „. dG(l,k, u) . „ sinh a , . dS(u.x) . „ 
[2 Z- +4 J=-\n coth M - l-2S(n,x)] V) ' w ]fc + 8 ^- G(k, I, fi) — — — — - x = 0. (4.4) 



ft dk 



Ok 

r 



dx 



It is easy to see that Eq. (4.4) defines a constant of motion effective Hamiltonian, H e s{x,k, fi,l) of the system and 
of the dynamical system. This constant can be called the from Eq.(4.4) we get 



H e g(x,k, = —2 — — — cos/c 



sinh 4 

4 \p coth/i - 1-2 S(u, x)] G(l,k,u). 
M 



(4.5) 



Again, it is easy to see from Eq.(4.4) that 



and k = — 



dk 



dx 



(4.6) 



together yield the original equations of x and k. In the 
subsequent analysis, the I = 1 case being the physically 
most relevant one, is analyzed. In another simplification, 
S(p, x) is approximated by its the most dominant first 
term. This simplification will not bring any qualitative 



change in the dynamics. 



B. The analysis of fixed point 

To obtain the fixed points of the set of equations, 
namely Eq.(4.2) and Eq.(4.3), we set x s = = k s . 
This then gives (x s ,k s ) = (|,ri7r) where p and n — 



0, ±1, ±2, 



To find the phase portraits around these 
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fixed points, we define z\ = x — x s and z 2 
We define for convenience and clarity 



We further define 
Bi(p,,gx) 



sinh n 
H 

Aq(ji) (p, coth^i — 1) 

4>(/*) — 

sinh^ 



(4.7) 



i4 x (/i) - (-If 2A 2 ( M ) 
(-l)"A)(/x) + 4 gi Bo(/i,p) 
16 ^ 2 9l A 2 (p). (4.8) 



Now, because of our stated assumptions, we get 

z\ = 2 Bi(fi,gi,p,n) z 2 
z 2 = (-If 2 B 2 (p, 9l ) zx. 



(4.9) 



This set of equations, in turn yields 



Bi(n,gi,p,n) 



(-If 



5 2 (M,3i) 



= Constant. (4.10) 



From Eq.(4.10), it is transparent that the system has only 
two types of fixed points. These are centres and saddle 
hyperbolic fixed points. We give below some in depth 
analysis of fixed points. 

We note that Ai(p), i = 0, 1, 2 are even fuctions of \i. 
This, in turn implies that Bi, i = 0, 1, 2 are also even 
functions of fi. We take gx > in this analysis. There 
will be no loss of generality by this assumption. We now 
note the followings. 

(i) B (a,p) — > 0, as p — > 0. It can also be shown 
that B^(p,p) is a monotonically increasing function of p. 
In other words, it is a monotonically increasing positive 
semidefinite function of p. These properties are true, ir- 
respective of p being even or odd. 

(ii) Bi(fi,gi,p,n) — > (—1)", as p. — ► 0. Furthermore, 
it can be seen by plotting this function against p that 
-Bi(/i, gi,p, n) > 0, when p — > oo. Both properties 
are true irrespective of the nature of p. Furthermore, 
when n is even, B\ (p, gx , p, n) is a monotonically increas- 
ing as well as a positive definite function of /i. On the 
other hand, when n is odd, there exists a value of /i, 
say fi r (gi,p) such that Bi (fj, r , g\ , p, n) = 0. Again, 
VriguP) nas the following properties, (a) For a given 
value of gi, p, r (jgi, evenp) > /U r (<7i, oddp), (b) irre- 
spective of the nature of p, the value of \i r monotonically 
decreases with increasing gx, and (c) p r (g\, evenp) — > 
/i r (gi, oddp) as g\ — * oo. See also Fig.l. 

(iii) B 2 {p,,gi) is a monotonically increasing positive 
semidefinite function of /i. 

We now consider case by case to determine the na- 
ture of fixed points from the linearized analysis. We 
note that the classifications are done using standard 
analysis pIL p2|. 




FIG. 1: This shows the variation of fi r { 
gi where fi r is the root of B\ (/i, gi , p, n) 
n is odd. I: p is odd and II: p is even. 



7i , p) as a function of 
Eq.(4.8) in the text. 



Case I : n and p are both even. In this case, both 
B\ (fi, g\ , p, n) and B 2 (fi,gi) are positive. Eq.(4.10) de- 
fines a hyperbola. So, the fixed point is a saddle point. 
Casell ; n is even, but p is odd. Since, the nature of 
these two functions remains the same, Eq.(4.10) defines 
an ellipse. So, we have a center. 

Case III: Both n and p are odd. Again, a > 
/i r (gx , evenp) . In this case, both B\ and B 2 are positive. 
So, Eq.(4.10) defines an ellipse and we have a center. 
Case IV: n is odd, but p is even. Furthermore, u > 
fj- r ( 9 i, evenp). In this case also, both B\ and B 2 are pos- 
itive. But due to the evenness of p, Eq.(4.10) defines a 
hyperbola and we have a saddle point. 
Case V: Both n and p are odd, but a < u r (gx, oddp). 
In this case, while B 2 is positive, Bx is negative. So, 
Eq.(4.10) defines a hyperbola and the fixed point is a 
saddle point. 

Case VI: n is odd, but p is even. Again, u < 
p, r (gx,oddp). In this case, Eq.(4.10) defines an ellipse. 
So, we have a center. 

Case VII: Both n and p are odd. Again, /j r (gx,oddp) < 
fi < fi r (gx , evenp) . In this case, both Bx and B 2 are 
positive. Since, p is odd, Eq.(4.10) defines an ellipse. So, 
we have a center. 

Case VIII: n is odd, but p is even. Furthermore, 
fi r (gx,oddp) < a < u r (gx, evenp). In this case, while 
£>i is negative, B 2 is, however, positive. Since, p is even, 
Eq.(4.10) defines an ellipse. So, we have a center. 
Our analysis is tabulated below. By Region I, it is meant 
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that /i < fi r (gi, oddp). Region II and III are respectively TABLE I : Nature of fixed points from linearization 
characterized by p r (gi , oddp) < fj, < fi r (9ii evenp) and analysis 
/i > n r (gi,evenp). 



r — £ 


fc s = n 7r 


Nature of fixed points 


P 


n 


Region I 


Region II 


Region III 


even 


even 


saddle 


saddle 


saddle 


odd 


even 


center 


center 


center 


odd 


odd 


saddle 


center 


center 


even 


odd 


center 


center 


saddle 



C. A numerical analysis of phase portraits 

Regarding the numerical investigation of phase dia- 
grams around fixed points, we note that there are all 
together six fundamental possibilities. Two possibilities 
arise from k s being even or odd. Again, for each case, 
there are three possibilities, depending on the magnitude 
of fi. In this discussion, we consider k s = and it. 
Along the line, k s — 0, we have a set of fixed points at 
x s = 0, ±i, ±1, ±|, ±2, • • •. According to the analy- 
sis, based on linearization x s — 0, ±1, ±2, • • • should be 
saddle hyperbolic fixed points. Between every two con- 
secutive saddle points, one should naturally expect cen- 
ters. This is perfectly borne out in the numerical inves- 
tigation (Fig. 2). It is also to be noted that Fig. 2 shows 
the phase diagram in the region I (/i < /i r ( <?i , oddp) ) . To 
understand the origin of a saddle fixed point between two 
consecutive centers, we note that two outermost ellipses 
encircling two centers will definitely touch at a point on 
the x-axis. Then, the flow of phase curves emanating 
from this point as well as in the vicinity of it will in- 
dicate a saddle hyperbolic fixed point. By the similar 
type of argument, the appearance of a center between 
two consecutive saddle fixed points can be understood. 
It is also found that the absolute critical value of k, below 
which localization of soliton occurs, increases monotoni- 
cally with increasing p. This behavior is consistent with 
the physics of the problem. 

We consider now the case of k s = it. It has certain 
interesting twists, which is quite a common occurence 
in the dynamics of nonlinear systems [p2[. In the region 
I, that is, pi < ^L r {gi,oddp), according to our stan- 
dard linearization based analysis, x s = 0, ±1, ±2, • • • 
should be centers. On the other hand, x s — ±i, ±|,--- 
should be saddle hyperbolic fixed points. The reason to 
have a saddle hyperbolic fixed point between two consec- 
utive centers has already been put forward in the pre- 
vious paragraph. Our numerical investigation also con- 
firms this result (Fig. 2). In the region II (p r (gi, oddp) < 
pi < p r {gi, evenp)), the linearization analysis tells us 
that all fixed points along the line k s = tt are centers. 
To understand this scenario, we note the following. Con- 
sider two centers, which are next but one to each other. 
Outermost ellipses, encircling these fixed points can also 
intersect. We then have a convex lens type region in the 
phase space (Fig. 3). Within this region of phase space, 
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4.00 




-1.00 ' 1 1 1 1 1 1 1 1 1 — 1 

-1.00 -0.50 0.00 w 0.50 1.00 1.50 



FIG. 2: The bottom part of the figure is the phase portrait 
that shows the trapped and the moving solitons in Region 
I for the fixed points having k 3 = 0.0. For moving solitons 
kmo = and ^- respectively. The upper part of the figure 
is the phase portrait that shows the trapped and the moving 
solitons in Region I, but for the fixed points having k s = ir. 
For the moving soliton k m7I — {x s } are shown in the 

figure. I — l,gi = 0.5, and fi — 1.0 for all curves in the 
figure. Eq.(4.5) is used to obtain this phase portrait. See also 
the text. 




2.50 1 

-0.70 -0.10 0.50 v 1.10 1.70 



4.20 




2.20 

-0.50 0.00 0.50 X 1.00 1.50 

FIG. 4: This phase portrait shows the trapped as well as the 
moving solitons in Region III for fixed points having k s — n. 
{x s } are, of course, shown in the figure. As per the text 
I = f. Furthermore, gi = 0.5, and fi = 1.2. For the 
moving soliton k mTr = ^ and ^p. 



phase curves will have no other option but to close around 
a point on the x-axis. This fixed point will also then be 
another center. This is clearly seen in Fig. 3, which de- 
scribes the phase diagram of the soliton dynamics in the 
region II. 

It is needless to say that the phase diagram is very 
interesting and unique. To the best of my knowledge, 
this is the first case where a chain of centers with no 
hyperbolic fixed point is found. The phase diagram in 
the region III is shown in Fig. 4. 

This region is defined by /i > ji r {g\,evenp). Accord- 
ing to the linearization analysis, along the line k s = n, 
x s = ±|, ±|,--- should be centers. This is indeed 
found in the numerical investigation. Again, the same 
linearization analysis tells that x s = 0, ±1, ±2, • • • will 
be saddle hyperbolic fixed points. From the numerical 
study with fj, = 1.2 and g\ = 0.5, we find that upto 
a certain distance from each of these fixed points, phase 
curves tend to diverge from the fixed points, indicating 
their hyperbolic structure. However, after a critical value 
of x, phase curves turn back to make close curves. This 
happens upto a certain maximum value of x from the rel- 
evant fixed point on the line k s = ir. After that value of 
x, phase curves close around the nearest center (Fig. 4). 
In case of x s = 0, it is found that the phase curve 
with the intial value of x, xq — 0.306 closes around 



FIG. 3: This phase portrait shows the trapped and the mov- 
ing solitons in Region II for fixed points having k a — n. {x 3 } 
are, of course, shown in the figure. As per the text I — 1. 

Other parameters, namely g\ = 0.5, and /i = 1. 15. For the 
„,„,,;„„ „„l;j.„„ u — 13vr 
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x s = 0. But, for xq = 0.307, the phase curve encircles 
x s = 0.5. So, for these two fixed points, the transition 
point is somewhere in between these two values. 

It is important to realize that in this part of the prob- 
lem we find three important phenomena that are encoun- 
tered in nonlinear dynamics, namely (i) bifurcation, (ii) 
discrete breathers and (iii) relaxation oscillations) 51, 52| , 

Hr{gi,oddp), that is, 



|53| , p4j . We note that for /i > n r {gi, oddp), that is, m 
regions II and III, all saddle points are changed to cen- 
ters (Fig. 3 and Fig. 4). So, we obtain a bifurcation in 
the phase portrait at ji = fi r (gi,oddp). Furthermore 
it is a saddle-center bifurcation. Again, at this bifur- 
cation point, it is easy to show that we have spatially 
localized but time periodic solutions. These localized so- 



lutions can be pinned at one of the permissible values of 
x s , namely at x s = ±i, ±|, • • •. These are by defini- 
tion discrete breathers [Q. To see it further, we examine 
Eq.(4.9). Since at /i = ^(gi, oddp), Bi = 0, from 
the first equation of Eq.(4.9) we find that Z\ — 0. The 
consistent solution is then z\ = 0. This in turn gives 
x = x s (oddp) = (h + |) where Zi = 0, ±1, ±2 
Then the second equation of Eq.(4.9) gives z\ = 0. 
This implies that k = k s — (2I2 + 1) ir where 
1-2 = 0, ±1, ±2, • • •. For this particular case , k s = it 
or I2 = 0. From the perturbative calculation of a ap- 
pearing in Eq.(3.1)|2^], it can be easily shown that a is 
proportional to r = Jt. The energy of the breather in 
terms original variables [Ml is 



J 



E[fi r (gi,oddp)} 



E, 



breather 



smh[i r ( g 1, oddp) 
1 + 2(7! 7 rrr — [1 — fi r (gi,oddp) coth /i r ( g\ , oddp) 



H[Hr{gi,oddp)} 



fi r (gi,oddp) 
4J 



A 



sinh [i r (gi, oddp) E[(i r (gi, oddp)]. 



E <-')'- 



fi r {gi,oddp) 



s = 1 



sinh 



fi r (g 1 ,oddp) 

(4.11) 



r 



We also note that for j- > 0, this is the minimum 
energy spatially localized state. Inasmuch as Eq.(2.5) 
has a reflection symmetry, another discrete breather will 
be obtained by the transformation, g\ —* ~gi- This 
is discussed later. We again note that our linearization 
analysis suggests that at /i = ^(gi^evenp), we have 
another bifurcation. This time it is a center saddle bifur- 
cation. However, we find that in the region III the phase 
portraits around x s = 0, ±1, ±2 • • • are not that of 
saddle points. In stead these phase diagrams consisting 
of closed orbits, are made up of both stable and unstable 
manifolds [M. So, here we have a relaxation oscillations 



dynamics as seen in Van der Pol oscillators [^l| in the 
Zeeman models of heartbeat and nerve impulse[p6| . We 
further note that in the present problem any relaxation 
oscillations dynamics is not physically allowed. So, these 
parts of the phase space will remain unaccessible to the 
system. 



STATIONARY LOCALIZED STATES 



J 



To study stationary localized states of Eq.(2.5), we seek for an oscillatory solution in the form[pl| 
^m( T ) = Qm( T ) exp[i(fcm — lot + do)] cxp[— 2ir] , to € Z 
where Q m is real and 00 is a constant phase factor. From the real and imaginary parts of Eq.(2.5), we have 

(n<9) m = luQ„ 



(5.1) 



cosfc (1 + Q m ) (Q m+ i + Q m -x) 
I 



2 (1+ Q 2 m ) Qm 9jCOS 2 kj (Q 



2 

m+j 



Q 



m-jj 



= 0; 



(5.2) 



Q m + Sinfc (1 + Q 2 m ) {Q m +l - Qm-l) 

I 

= -2 (1+ Q 2 m ) Q m Y 9j cosk i sinfcj (Q 
3 = 1 



2 



i 1 X 
^ m—j i 



(5.3) 



r 



where Q is the column vector (Qi,Q%,'-' Q r , 



■) and f2 is the matrix defined by the left hand side of 
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Eq.(5.2). Eq.(5.2) with vanishing boundary condition 
constitutes a nonlinear eigenvalue problem for localized 
states^, [ll], [$l). Eq.(5.3) determines the time evo- 
lution of the localized states. It is a trite algebra to 
show that when gj = for all j, Q, m (r) given by 
Eq.(3.1) is a solution of Eq.(5.2) and Eq.(5.3). When 
k = or 7r, Q m , m £ 2 is stationary. A localized 
state is called "staggered" if k — w, and " unstaggered" 
if k = 0. We further observe that under the transforma- 
tion, Q m — ► (— l) m <5m, Eq. (5.2) will remain invariant if 



lj — > — u> and gj — > — gj j G I. This is indeed in ac- 
cordance with the reflection symmetry of Eq.(2.5). This 
result in turn implies that if an unstaggered localized 
state, Q m exp[— iu>r] is a solution of Eq.(5.2), the corre- 
sponding staggered localized state, (— l) m Q m exp[i(x>r] is 
then a solution of the same, provided gj, j S / changes 
sign for all j. To determine locations of these localized 
states, we shall specialize on I = 1 case. This simplifi- 
cation will not affect the merit of the discussion. From 
Eq.(5.2), we have 



J 



-2 cos A; — (cosfc + 2g\ cos 2 k) 



Era Qm 



— 2gi cos 2 k 



Era Qm (Qr 



Em Q r 



Qm-l) 



(5.4) 



We suppose that Q m > 0, m e Z and \gi\ <C 1. 
Then, the staggered state lies above the phonon band, 
while the unstaggered state lies below. We further note 
that when A — > 0, \g\\ — > oo. In this case, there is 
no localized state, staggered or unstaggered, below the 
phonon band if g\ < 0, or above the phonon band for 
.91 >0. 

The next important point is the effective mass of un- 



staggered and staggered localized states^!], |3l). For this 
we consider the Hamiltonian, H c ff given by Eq.(4.5). This 
consideration will allow us to obtain expressions for the 
effective mass for almost unstaggered and almost stag- 
gered states along with the fully unstaggered and fully 
staggered localized states [21|. To calculate the effective 
mass of nearly unstaggered localized states, we let k — > 
in Eq.(4.5). This in turn yields 



J 



lim H e s(x,k,/j,) 

k — > 



-%{x) k 2 - 2A (p)- AgMn) + 8 9l A 2 {p) cos2nx + 0(fc 4 ), 



(5.5) 



r 



where we define 

m~g(x) = 2A Q (n) + 8g 1 A 1 (fi) - 16 g 1 A 2 (n) cos27ra; 

(5.6) 

and Ao, Ai, and A2 are already defined in the text (see 
Eq.(4.7)). We note that m c s(x) is even function of both 
fi and x. It is to be noted that because of the pro- 
posed spatial dependence of m e fi, this is a generalized 
definition of the quantity. It can be easily seen that for 
gi positive semidefinite, m e s(x) is positive definite for 
/i G [—00, 00]. So, when the coupling constant, g\ > 0, 



the effective mass of unstaggered and nearly unstaggered 
states is positive definite. 

To obtain the effective mass of staggered localized 
states, we put k = it — 9 and then let 9^0. 
This procedure gives 

-m~g(x) = 2 A (fx) - 851 A±(ji) - 16 51 A 2 (/x) cos27ra. 
Consequently, we get from Eq.(5.7) that when fi — > 



0, 



-2. On the other hand, in the limit 



jj, ^> 1, the asymptotic form of Eq.(5.7) is 



J 



^EM [l-2 9l expM - 2 5l ^M(l + 2cos27rx)]. 



(5.8) 



r 



Let us assume that 2 g\ — exp[— (1 — e)/i], e > and 
e n — > 0. We then see from Eq. (5.8) that m e g(x) changes 
sign. So, depending on the value of g%, there will be a 



critical value of /1, [i CI (gi), such that for /1 > fi CT (gi), 
staggered as well as nearly staggered localized states will 
have positive effective mass. For a more precise analy- 



11 



4.0 




0.0 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.4 0.8 Qi 1.2 1.6 2.0 

FIG. 5: This shows the variation of /j} CI and /u£ r as a function 
of gi . Curves I and II respectively. These critical points deter- 
mine the signature of the effective mass of staggered localized 
states. Eq(5.7) in the text. 



sis of the effective mass, we note that Eq.(5.7) has two 
critical values of fi, namely pr CI , and /x" r , depending on 
whether cos 2irx = 1 or — 1. Again, these two critical 
values of \i move to each other as g\ increases. This re- 
sult is indeed substantiated by the numerical analysis of 
Eq.(5.7) as shown in Fig. 5. It is again to note that for 
/i > /i" r , the effective mass is positive definite, while for 
\i < /ij. r , the effective mass is negative definite. This in 
turn implies that the system will have a strongly localized 
staggered state for > iMcrCSi)!- We further note from 
Eq.(3.1) that ijT 1 gives the measure of the localization 
length. So, states for which \i — ► in Eq.(3.1) have long 
localization lengths. We conclude from this analysis then 
that stationary staggered localized states in the vicinity 
of upper phonon band edge will have negative effective 
mass. 



VI. A NUMERICAL STUDY OF EQ.(2.5) 

A. A study of soliary wave like solutions 

For the purpose of numerical integration, we first re- 
place \& m , m = 1,2,3, • • • by * m exp(-2ir) in Eq.(2.5). 
Next we write * m = *i m + 3 ^ 2 m- From Eq.(2.5), 
we then get 



*lm - (1 + *L + *2m)[(*2(m+l) + *2(m-l)) + 2 ^ 9j { (* lm * 1(m+J -) 

3=1 

+ ^2m^2(m+j))^2(m+j) + (*lm*l(m-j) + *2m*2(m-j))*2(m-j)}] 

I 

*2m = (1 + *L + *lm)[(*l(m+l) + *l(m-l)) + 2 £ ft {(*l m *l( m+J ) 

3=1 

+ *2m V l'2(m+j))*l(m+j) + (*lTn*l(m-j) + *2m*2(m-3))*l(m- 3 )}]- 



(6.1) 



(6.2) 



r 



For the initial condition, we use[|| 

sinh u 

*lm( T = 0) = ; cos km 



^ 2m (r = 0) 



COsllTO^ 

sinh fx 
cosh m/x 



sin fern. 



(6.3) 
(6.4) 



The Fourth order Runge-Kutta method is used to inte- 
grate these equations. The chosen time interval for all 
calculations is 10~ 4 . Furthermore, Eq.(2.4) with A = 1 
is used to check the accuracy of the integration. In this 
study an one dimensional lattice, comprising of 257 lat- 
tice points with unit lattice spacing is used. The initial 
condition is centered about the middle of the lattice, that 
is about the 129th site. For the numerical analysis, we 



consider two cases, namely (i) I = 1 and (ii) I = 2 but 
gi — in Eq.(2.5). For all cases that are presented in 
this section \x = 1 and gj — 0.5, j — 1, 2. 

We note that in the first case the permissible values of k 
are ±5 while for the second these values are ±-| and ±^r. 
For the first case the absolute amplitude of the solitary 
wave like solution (Q m , see Eq.(3.1)) as a function of 
site and time (m, r) is shown in Fig. 6 for k = -|. 
Furthermore, for this case the time evolution of the ab- 
solute amplitude of the solitary wave (Qto(t)) for some 
arbitarily chosen sites, m are shown in curves I, II, III in 
Fig. 7. For the second case k = — ^ is chosen and the 
space-time evolution of the Q m (r) is shown in Fig. 8. We 
note that in all cases the initial profile moves undistorted. 
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FIG. 6: The N-AL dynamics of an initial Ablowitz-Ladik soli- 
ton with k = tj, Eq.(2.5). Furthermore, I = 1, gi = 
0.5 and /i = 1.0. The number of sites in the chain is 257 and 
the origin is taken at the center of the chain. 



FIG. 8: The N-AL dynamics of an initial Ablowitz-Ladik soli- 
ton with k = — ^f, Eq.(2.5). Furthermore, 1 = 2, gi = 
0.0, Q2 = 0.5 and fi = 1.0. The number of sites in the 
chain is 257 and the origin is taken at the center of the chain. 




0.0 2.0 4.0 x 6.0 8.0 10.0 

FIG. 7: The N-AL equation (Eq.(2.5)) determined time evo- 
lution of Q m (i~) (Eq.(3.1)) for some chosen values of sites, m 
. k = and ^ respectively. For both cases 1 = 1, gi = 
0.5 and /i = 1.0. The number of sites in the chain is 257 and 
the origin is taken at the center of the chain. I, 1(a): m = 
130, II, 11(a): m = 140, III, 111(a): m = 150. The (a) series 
is for k — ^ while the other one is for k = ^. 



These pictures are expectedly identical to the space-time 
evolution of the pure A-L solitons. This is also verified 
by numerical integration. 



B. A study of the PN pinning of solitary waves 

Another important aspect that we study numerically 
here is the effect of PN potential, which arises due to 
the nonintegrability term in Eq.(2.5) on the space-time 
evolution of the initial profile given by Eq.(6.3) and 
Eq.(6.4)||l|, H). We note that the nonintegrability effect 
is maximum when k = in Eq.(3.1). On the other 
hand it totally vanishes at |fc| = ^ for I — 1. We study 
again the propagation of an initial profile for I = 1, 
but with /i = 1.0 and gi = 0.5. The accuracy of the 
integration is checked by the constancy of A/"(Eq.(2.4)). 
It is found that the loss of constancy becomes more 
and more discernible as \\k\ — §| —* §• This makes 
the simulation for small values of k less reliable. Fig. 9 
shows the space-time evolution of Q m (r) (Eq.(3.1)) for 
k = J. We note that the dispersion and the distortion 
of the initial profile is very transparent in the figure 
(Fig. 9). For better understanding, we also show through 
curves 1(a) to 111(a) in Fig.7 and 1(a) to IV(a) in Fig. 10 
the time evolution of the Q m (r) for some arbitrarily 
chosen site or m values for k = | and j respectively. 
We note that the distortion of the initial profile and the 
intensity of the phonon tail increase as k decreases from 



13 



Qm(t) 




| to |. This is, of course in the expected direction. 
Furthermore, in both cases the initial profile slows down 
by leaving phonon behind in the propagation. Through 
curves I to IV in Fig. 10 the time evolution of Q m (r) 
for same values of m from the Ablowitz-Ladik equation 
are included for comparison. This slowing down is very 
transparent in the figure (Fig. 10). So, it is reasonable 
to assume that these moving profiles will ultimately be 
trapped. However, we do not observe the total trapping 
in our simulations. We possibly need large integrability 
breaking term to see the total trapping in the time of 
simulations. In case of A: = ? , the nonintegrability term 
in Eq.(2.5) vanishes, and expectedly the initial profile 
(Eq.(6.3) and Eq.(6.4)) propagates undistorted with- 
out any change in its speed and leaving no phonon tail 
behind as can be seen in Fig. 6 and curves I to III in Fig. 7. 



FIG. 9: The N-AL dynamics of an initial Ablowitz-Ladik soli- 
ton with k — -|, Eq.(2.5). Furthermore, I = 1, gi = 
0.5 and /i = 1.0. The number of sites in the chain is 257 and 
the origin is taken at the center of the chain. 



1.50 




0.00 2.00 4.00 x 6.00 8.00 10.00 

FIG. 10: Same as Figure 7, but with k — -| for both cases. 
I, 1(a): m = 130, II, 11(a): m = 135, III, 111(a): m = 140, IV, 
IV(a): m = 142. While (a) series is for the N-AL equation, 
the other one is for the Ablowitz-Ladik soliton. 



C. A study of interaction of two solitary waves 



An important and also extensively studied field in non- 
linear dynamics is the investigation and the understand- 
ing of interaction of two or more solitary waves of nonin- 
tegrable nonlinear equations and the components of vec- 
tor solitons in integrable nonlinear equations. So far 
only continuous nonlinear nonintcgrable equations are 
studied^, |§ §§, |§ §(| || || g§. This is due 
to the nonavailability of any discrete nonlinear equation 
having solitary wave like solutions. In this context, there- 
fore, the N-AL equation (Eq.(2.5)) that is proposed here 
assumes a very great significance. This equation gives 
us the opportunity for the first time to study the inter- 
action of solitary waves of a discrete nonlinear equation. 
To investigate th is p roblem, as the initial condition we 
take the function^, |3(| 



14 



T , n \ ri sum u , 

cosn[/i(m — xq)J 



l^mCOl 




FIG. 11: The N-AL collision dynamics of two initial A-L 
pulses, = — fci = f . Z = 1, 3i = 0.5 and /x = 1.0 
as mentioned in the text. The number of sites in the chain is 
313 and the origin is taken at the middle of the chain. 

We note that the function has two peaks at x = Xq and 
x = — xq respectively. The velocity of the peak at xq 
is sm ^ 1 M sin ki while the velocity of the other peak at 

x = —xq is am sin k 2 ■ We further note that when 
|a?o| ~~ ¥ °°i Eq.(6.5) gives two A-L type solitons 
(Eq.(3.1)). This analysis is also done numerically using 
the fourth order Runge-Kutta method, and the origin is 
placed at the center of the one dimensional chain. For all 
studies, we take 2|xo| = 30 units, and /i = 7*2 = 1-0 
The first one is as usual the I = 1 case. For this anal- 
ysis we take for all cases /i = 1.0 g\ — 0.5. Inasmuch 
as \k\ — S in Eq.(3.1) is an exact one soliton solution of 
Eq. (2.5) , we consider first the case k 2 — — fci = 5 . In 
this case, we find that two peaks emerge after the colli- 
sion without any change in shape and without any emis- 
sion of phonon (see also ref.34). Of course, there can be 
phase shifts in these peaks after collision, but this is not 
investigated. This case of collision of solitons is shown 
in Fig.ll. We further note that the result is found to be 
independent of the value of g%. But, we did not check 
the effect of changing fi on the collision. The other case 
that is considered is fci = — jfjg and k 2 = § ■ When 
9i = 0, we find that two solitons collide and after the 




FIG. 12: Same as in Fig.ll, but with ki = and ki = 
— 5 . I — 1, gi — 0.5, and /x = 1.0. In this simulation, the 
number of sites is 257. 



collision again two peaks emerge without any change in 
shape and without any emission of phonon. This is not 
shown here. On the other hand, when we give a non-zero 
value of gi (in our case 0.5), we find the fusion of two 
solitons on collision. Inasmuch as k 2 = the direction 
of the maximum velocity is in the direction of increas- 
ing lattice sites. The fused solution expectedly moves in 
that direction. But, at the same time it emits phonons, 
as shown in Fig. 12 

The second case that we study is the case of I — 2, but 
gi = 0.0. In this case we have four permissible values 
of k, namely k — ±^p, and ± j. Of course, there are 
only two velocities, just like the previous one. However, 
if choose k\ and k 2 from these allowed values of k, for 
\xq\ — > oo, Eq.(6.5) will be the solution of Eq.(2.5). In 
this analysis, we choose [i = 1.5, g 2 = 0.5, and2|a;o| = 
30 units. For this case the space-time evolution of a soli- 
tary wave with k = — ^ is already shown in Fig. 8. 
In the first case we take k\ = —k 2 = — ^ and 
the numerical simulation is shown in Fig. 13. In this 
cae we see that two solitary waves do not pass each 
other. On the other hand after coming close to each 
other, they are repelled. Other cases that are studied are 
fci = - f , but k 2 = ^, andfc! = - 2f, but k 2 = \ 



15 




190 



FIG. 13: The N-AL collision dynamics of two initial A-L 
pulses. k 2 = -hi = ^f- I = 2, gi = 0.0, g 2 = 
0.5 and n = 1.5 as mentioned in the text. The number of 
sites in the chain is 313 and the origin is taken at the middle 
of the chain. 

respectively. Our numerical simulations for these cases 
are shown in Fig. 14 and Fig. 15 respectively. In the first 
case, we see that two solitary waves fuse after collision 
and then move in the direction of positive velocity si- 
multaneously emiting phonons. In the second case, how- 
ever, we see the total destruction of the initial profile, 
Eq.(6.5) after the collision. This implies that there is 
no definite pattern in the collision process. Furthermore, 
the result of collision appears to be very sensitive to the 
phases of the colliding solitary waves. So, there should 
be attempts to understand this collision physics analyt- 
ically, using perturbative method and collective coordi- 
nate method 

VII. SUMMARY 

A new class of nonintegrable Hamiltonians with tun- 
able nonlinearities is propsed here to derive a class of 
(1 + 1) dimensional nonintegrable discrete nonlinear 
Schrodinger equations, which include both some known 
nonlinear equations such as Ablowitz-Ladik nonlinear 
Schrodinger equation §, §, IN-DNLSg, ||, [«§ and a 
new subset, collectively christened as the N-AL equa- 
tion(Eq.(2.5)). The relevant equations are obtained from 
the proposed Hamiltonian by the standard procedure of 
classical mechanics, but by employing a standard gen- 
eralized definition of Poisson brackets, as shown in the 
Appendix A below. In this paper two cases of the N-AL 




FIG. 14: Same as in Fig. 13, but with ki = — f , and 
k 2 = . I = 2, gi = 0.0, g 2 = 0.5, and fj, = 1.5. In 
this simulation, the number of sites is 313. 




FIG. 15: Same as in Fig. 14, but with ki = — '^f, and = 
j. I = 2, 31 = 0.0, g 2 = 0.5, and fi = 1.5. In this 
simulation, the number of sites is 313. 
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equation are investigated analytically as well as numeri- 
cally. These cases are (i) I — 1, and (ii) I — 2, but gi = 
0.0 in Eq.(2.5). 

An important characteristic of the N-AL equation is 
shown to be that it can allow spatially localized states 
of A-L type (Eq.(3.1)) for certain permissible values of 
the parameter k to travel without distortion and without 
cmiting phonons. This is found both analytically as well 
as numerically. So, these special solutions have in them 
the required balance of the nonlinearity and dispersion]^] 
and also are transparent to the PN potential arising from 
the lattice discreteness [|$l| fil) . The transparency of these 
solutions to the PN potential proves that these have the 
continuous translational symmetry of the one-soliton so- 
lution of the Ablowitz-Ladik equation (Eq.(3.1))|pl|. 

Trapped and moving solitons are found from Eq.(2.5). 
This analysis is done using a standard perturbative 
procedure pl|, |22]] . The most interesting part of this anal- 
ysis, however, is that it shows the existence of saddle 
center bifurcations. This in turn implies the presence of 
minimum energy breathers in the system. Breathers are 
by definition time periodic, spatially localized solutions 
of equations of motion for classical degrees of freedom 
interacting on a lattice jH^] . We note that eq.(4.1) satisfy 
these requirements. Since, the saddle center bifurcation 
is seen in the perturbative method, it is, therefore, im- 
perative to study the existence of breathers in the system 
numerically. We further note that in this system we also 
find a chain of centers. The origin of this situation is ex- 
plained in the text. One interesting consequence of this 
occurence is phase portraits showing relaxation oscilla- 
tions in regions, which are physically forbidden parts of 
the phase space for this system. Since phase portraits 
are very rich in characteristics, it is worthwhile to study 
them quite elaborately by varying extensively parame- 
ters, such as jji and g\ and and also including more terms 
in S(n,x) (Eq.(4.1). 

The presence of stable localized states in Eq.(2.5) is 
also investigated, albeit a more systematic analysis, pre- 
sumably by the discrete variational method is required. 
This work is indeed in progress. The effective mass of lo- 
calized states are analyzed here using the effective Hamil- 
tonian, given by Eq.(4.5)pl|. It is found that nearly 
unstaggered as well as unstaggered localized states will 
always have positive effective mass. On the other hand, 
nearly staggered and staggered localized states can have 
both negative and positive effective depending on the 
value of n in Eq.(3.1). However, large width but small 
amplitude nearly staggered and staggered localized states 
do have the expected negative effective mass. 

It is further shown here numerically that all other soli- 
tary wave profiles (Eq.(6.3) and Eq.(6.4)) under the dy- 
namics described by Eq.(2.5) and for which the nonin- 
tegrable term in Eq.(2.5) does not vanish, suffer distor- 
tion and emit phonons in the propagation. Further in- 
vestigations along this line by altering the strength of 
the nonintegrable term by altering the value of k reveal 
that the extent of distortion and the extent of emission 



of phonons by the initial solitary wave profiles (Eq.(6.3) 
and Eq.(6.4)) depend on the magnitude of the noninte- 
grable term. We conclude from these two findings that 
any general solitary wave profile due to the presence of 
the nonintegrable term in Eq.(2.5) is subjected to both 
the dynamical imbalance between the nonlinearity and 
the dispersion, and to the PN potential. We note that 
the continuous translational symmetry of these solitary 
wave profiles breaks down due to the nonintegrable term 
in Eq.(2.5). Consequently, these profiles scan the lattice 
discreteness in the propagation. This, in turn implies 
that the moving profiles experience the effect of PN po- 
tential in the dynamics. Since, any general solitary wave 
profile suffers from these two factors, we expect that the 
profile can either be trapped by emiting phonons or the 
localized structure will spread uniformly due to the ef- 
fect of extra dispersion. Notwithstanding the slowing 
down of initial profiles are observed in our simulations, 
our numerical study in this matter, however, is not con- 
clusive. So, a more elaborate study increasing the cou- 
pling constants in the nonintegrable term in Eq.(2.5) and 
also changing the parameter /i in Eq.(3.1) is needed. This 
will be done in future. 



Finally, the N-AL equation (Eq.(2.5))is the discrete 
analog of the famous c/> 4 equation, however with an im- 
portant difference, that this equation only gives dynam- 
ical solitary wave like solutions. None the less, the im- 
portance of the N-AL equation in this regard cannot be 
overlooked. This equation gives the opportunity for the 
first time to study the collision dynamics of solitary wave 
profiles using a discrete nonlinear nonintegrable equation. 
Here only the collision of two solitary wave profiles are 
studied (Eq.(6.5)). In this study we find that at least 
the interaction of two solitary wave profiles have some 
universal features, meaning that the outcome does not 
depend critically on the nature of the nonlinear equa- 
tion. It is also equally true that we do not find so far 
any systematic pattern in the dynamics. It is, therefore, 
necessary to study this problem in details both analyt- 
ically and numerically to discern if there is any specific 
pattern in the dynamics. Another imporatnt aspect to 
study is the fractal structure of the emeging profile after 
the collision as a function of the initial velocity and the 
strength of the nonintegrable term j$7. 
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We conclude by noting that this study offers a very sig- 
nificant insight into the transport properties of the DHM. 
Furthermore, the N-AL equation can be used to study 
transport properties of localized states in soft molecular 
chains, having the coexsistence of nonlinearity and disor- 
der, nolinearity and quasiperiodicity or incommesurabily 
in various parameters, such as the hopping integral, J and 
the nonintegrable coupling constants , <7j , j — 1, 2, 
Furthermore, the aplicability of the N-AL equation in the 
study exciton dynamics in photosynthesis and in molec- 
ular solids should be explored. 
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APPENDIX A: GENERALIZED POISSON 
BRACKET 

We consider a dynamical system having 2 N general- 
ized coordinates, {(j> n , </>^}, n = 1, ...TV. Let U and 
V be any two general dynamical variables of the system. 
We now define a nonstandard Poisson bracket Q 

, dU dV dV dU w , , , , 
{U, V} { ,, r} = £ (or 7^ ~ 5CT TUT) (! + A W 



iV 



= £ (1 + A |0„| 2 ) ( 

n = 1 

= {U, H} { ; 



<9U <9H 



d(j)* n d<j) n 
(A6) 



d< Pr> 

When [/ = m and V = we then have [21. 31. ."> 

{4>m, </>*}{0, 0*} = (1 + A |<Am| 2 ) <W (A2) 
From Eq.(Al) we further have for {m, 1} 

{<t>m, M{4>, r} = {^m> 0*} = °- ( A3 ) 

Then when U — <p m and V = H, we have from Eq.(Al) 

<9H 



{</> m ,H} { i = (1 + \\<p m \ ) 



(A4) 



We then write for the dynamical evolution of the m-th 
generalized coordinate, <ft m fpl|, 



<9H 



dt 



{</> m ,H} w , j,*} = (1 + A |</> m | ) ^— . (A5) 



This is consistent because 

N 



dU 
~dT 



n — 1 T 



dU 



We note now that when U = H, we have 4^ = 0. In 

5 at 

other words, H is a constant of motion. Consider J\f given 
■by Eq.(2.4). To show that it is a constant of motion, we 
writer! = H + H a - §(£ - J) Af, where we define 



Hn 



Hi 



J E {4>n <]>n+l 
1 ' 

- o E E 9i 



<j>* 0„) 2 .(A8) 



n j = l 



We then find that 



l ir = A E^«7^- 



5H ffEi 



(A9) 

By using (A7) and ( A8) , it is very simple to show that the 
right hand side of (A9) is zero. Hence, Af is a constant 
of motion. 
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